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ABSTRACT 

We present detailed comparisons of the intracluster medium (ICM) in cosmological Eulerian cluster simula- 
tions with deep Chandra observations of nearby relaxed clusters. To assess the impact of galaxy formation, we 
compare two sets of simulations, one performed in the non-radiative regime and another with radiative cooling 
and several physical processes critical to various aspects of galaxy formation: star formation, metal enrichment 
and stellar feedback. We show that the observed ICM properties outside cluster cores are well-reproduced in 
the simulations that include cooling and star formation, while the non-radiative simulations predict an overall 
shape of the ICM profiles inconsistent with observations. In particular, we find that the ICM entropy in our 
runs with cooling is enhanced to the observed levels at radii as large as half of the virial radius. We also find 
that outside cluster cores entropy scaling with the mean ICM temperature in both simulations and Chandra 
observations is consistent with being self-similar within current error bars. We find that the pressure profiles of 
simulated clusters are also close to self-similar and exhibit little cluster-to-cluster scatter. We provide analytic 
fitting formulae for the pressure profiles of the simulated and observed clusters. The X-ray observable-mass 
relations for our simulated sample agree with the Chandra measurements to ss 10% - 20% in normalization. 
We show that this systematic difference could be caused by the subsonic gas motions, unaccounted for in X-ray 
hydrostatic mass estimates. The much improved agreement of simulations and observations in the ICM pro- 
files and scaling relations is encouraging and the existence of tight relations of X-ray observables, such as Yx, 
and total cluster mass and the simple redshift evolution of these relations hold promise for the use of clusters 
as cosmological probes. However, the disagreement between the predicted and observed fractions of cluster 
baryons in stars remains a major puzzle. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 



1. INTRODUCTION 

Clusters of galaxies are fascinating astrophysical objects 
and laboratories for studying galaxy formation and structure 
formation in general. At the same time, clusters can provide 
cosmological constraints that are complementary to those ob- 
tained with other methods such as temperature anisotropics 
of the cosmic microwave background, Type la supernovae, 
and weak lensing (e.g., Voit 2005; Tozzi 2006; Borgani 2006; 
Albrecht et al. 2006). Cosmological applications of clus- 
ters include cluster counts and their evolution with redshift 
(e.g., Henry & Arnaud 1991; Markevitch 1998; Ikebe et al. 
2002; Vikhlinin et al. 2003), spatial distribution (e.g., Miller 
et al. 2002), and the angular-diameter distance measurements 
(Allen et al. 2004; LaRoque et al. 2006). Detailed ob- 
servations of merging clusters provide unique insights into 
the physics of the intracluster plasma (e.g., Vikhlinin et al. 
2001a,b; Markevitch et al. 2003) and provide key evidence 
for the existence and properties of dark matter (Markevitch 
et al. 2004; Clowe et al. 2006). 

All cosmological applications of clusters, at least to a cer- 
tain degree, rely on solid understanding of the physics of 
their formation. Given that clusters are nonlinear collapsed 
systems, numerical cosmological simulations are the method 
of choice for their theoretical studies. Modern cosmologi- 
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cal codes using A/-body and numerical hydrodynamics tech- 
niques can accurately follow dynamics of dark matter and 
gaseous baryonic components in their full complexity dur- 
ing the hierarchical build-up of structures. Yet, more realistic 
modeling of clusters requires inclusion of additional baryonic 
processes. For example, to model formation of cluster galax- 
ies, we need, at the very least, to correctly treat energy dissi- 
pation due to radiative losses by baryons, and conversion of 
gas into stars. In addition, any feedback in the form of energy 
injection and metal enrichment from supernova winds (e.g., 
Metzler & Evrard 1994; Valdarnini 2003) and active galactic 
nuclei (Briiggen et al. 2005; Sijacki & Springel 2006; Cat- 
taneo & Teyssier 2007), and injection of non-thermal cos- 
mic rays at large-scale shocks accompanying cluster forma- 
tion (Pfrommer et al. 2007) can alter the thermodynamics of 
the intracluster gas. 

Although our understanding of details and relative impor- 
tance of these processes is currently sketchy, the simulations 
with specific assumptions about them are highly predictive, 
which should make models falsifiable. In particular, by com- 
paring observed cluster properties with the results of simula- 
tions that include various physical processes described above 
we can learn a great deal about these processes and their role 
in cluster formation. 

Over the last two decades, such comparisons were used ex- 
tensively to put constraints on the deviations of ICM ther- 
modynamics from the simple self-similar behavior, described 
originally by Kaiser (1986, 1991). The first studies of the ob- 
served correlation of cluster X-ray luminosity, Lx, and spec- 
tral temperature, Tx, unambiguously showed that its slope 
is steeper than the slope predicted by the self-similar model 
(e.g., Edge & Stewart 1991; Henry & Arnaud 1991; White 
et al. 1997; Markevitch 1998; Allen & Fabian 1998; Arnaud 
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& Evrard 1999). In addition, the slope of the Lx - 7x relation 
steepens for the lowest mass clusters (e.g., Helsdon & Pon- 
man 2000; Finoguenov et al. 2002, 2007). Deviations from 
self-similarity were shown to be the strongest in the cores of 
clusters (e.g., Markevitch 1998; De Grandi & Molendi 2002; 
Vikhlinin et al. 2006) and were widely interpreted as evidence 
for preheating of the intracluster gas by energy from super- 
novae and AGN feedback (e.g., David et al. 1991; Kaiser 
1991; Evrard & Henry 1991; White 1991; Wu et al. 2000; 
Bialek et al. 2001; Borgani et al. 2001, 2002; Nath & Roy- 
chowdhury 2002). 

Alternative explanation was proposed by Bryan (2000), 
who argued that cooling and condensation of the gas accom- 
panying formation of cluster galaxies can reduce the ICM gas 
density and increase its temperature and entropy to the ob- 
served levels (see also Voit & Bryan 2001; Voit et al. 2002). 
This explanation was borne out by cosmological simulations 
(Pearce et al. 2000; Muanwong et al. 2001; Valdarnini 2002; 
Dave et al. 2002; Kay et al. 2004, 2007). However, the amount 
of gas that condenses out of the hot ICM in cosmological sim- 
ulations due to cooling (e.g., Suginohara & Ostriker 1998; 
Lewis et al. 2000; Pearce et al. 2000; Dave et al. 2002; Ettori 
et al. 2004) appears to be a factor of 2-3 too large compared 
to the observed stellar mass in clusters (Lin et al. 2003; Gon- 
zalez et al. 2007). Thus, the X-ray measurements appear to 
be consistent with a large fraction of cooling gas, while the 
optical estimates of stellar mass indicate that this fraction is 
small. 

Modern X-ray observations with Chandra and XMM- 
Newton allow us to study the ICM properties with unprece- 
dented detail and accuracy. Their superb spatial resolution 
and sensitivity enable resolved, accurate X-ray brightness and 
temperature maps over a large fraction of the cluster virial 
radii. The X-ray measurements also enable accurate mass 
modeling of relaxed clusters with the assumption of hydro- 
static equilibrium of the ICM in the cluster potential. These 
observations can therefore be used for detailed comparisons 
of both global cluster properties and their profiles with simu- 
lation results, which provide more stringent tests for the mod- 
els of the ICM thermodynamics. In particular, such com- 
parisons can shed some light on the apparently contradictory 
lines of evidence as to the efficiency of cooling in clusters de- 
scribed above. 

To this end, in the present study we focus on the effects of 
radiative cooling and star formation on the observable prop- 
erties of clusters and compare results of simulations in both 
non-radiative and radiative regimes with the current X-ray 
data. Namely, we use two sets of simulations started from the 
same initial conditions. Both sets treat collisionless dynamics 
of dark matter and hydrodynamics of diffuse gas with high- 
resolution using the adaptive mesh refinement technique. In 
the baseline set of cluster simulations, the gas is modeled in 
non-radiative regime and thus does not reach high densities 
and is not allowed to form stars. The second set of simu- 
lations includes several processes accompanying galaxy for- 
mation: gas cooling, star formation, metal enrichment and 
thermal feedback due to the supernovae. Comparison of the 
simulated profiles in these two sets of simulations to those of 
observed clusters allows us to gauge the role of galaxy forma- 
tion in shaping properties of the ICM. As we show in § 4, the 
simulations that include galaxy formation processes provide 
a considerably better match to the observed ICM profiles out- 
side cluster cores compared to the non-radiative simulations. 

The paper is organized as follows. In § 2 we describe cos- 



mological cluster simulations. The methods used to analyze 
the simulations and the brief description of observations used 
in our comparisons are given in § 2 and 3, respectively. We 
present results of comparison of the ICM density, tempera- 
ture, entropy, and pressure profiles in § 4, and integrated quan- 
tities such as spectral X-ray gas temperature, gas mass, and 
pressure, in simulations and observations in § 5. We discuss 
our results and conclusions in § 6. 

2. COSMOLOGICAL CLUSTER SIMULATIONS 

In this study, we analyze high-resolution cosmological 
simulations of 16 cluster-sized systems in the flat ACDM 
model: Q m = 1 - Q A = 0.3, O b = 0.04286, h = 0.7, 
and <T8 = 0.9, where the Hubble constant is defined as 
100/! km s- 1 Mpc-',and an erg is the power spectrum normal- 
ization on an 8/z~' Mpc scale. The simulations were done with 
the Adaptive Refinement Tree (ART) A^-body+gasdynamics 
code (Kravtsov 1999; Kravtsov et al. 2002), an Eulerian code 
that uses adaptive refinement in space and time, and (non- 
adaptive) refinement in mass (Klypin et al. 2001) to reach the 
high dynamic range required to resolve cores of halos formed 
in self-consistent cosmological simulations. The same set 
of cluster simulations was used in our related recent studies 
(Kravtsov et al. 2006; Nagai et al. 2007), where additional de- 
tails can be found. We provide a description of the simulation 
details here for completeness. 

The Af-body+gasdynamics cluster simulations used in this 
analysis follow collisionless dynamics of dark matter and 
stars, gasdynamics and several physical processes critical to 
various aspects of galaxy formation: star formation, metal en- 
richment and thermal feedback due to Type II and Type la 
supernovae, self-consistent advection of metals, metallicity- 
dependent radiative cooling and UV heating due to cosmo- 
logical ionizing background (Haardt & Madau 1996). The 
cooling and heating rates take into account Compton heating 
and cooling of plasma, UV heating, and atomic and molec- 
ular cooling, and are tabulated for the temperature range 
10 2 < T < 10 9 K and a grid of metallicities, and UV inten- 
sities using the Cloudy code (ver. 96b4; Ferland et al. 1998). 
The Cloudy cooling and heating rates take into account metal- 
licity of the gas, which is calculated self-consistently in the 
simulation, so that the local cooling rates depend on the local 
metallicity of the gas. Star formation in these simulations was 
done using the observationally-motivated recipe (e.g., Kenni- 
cutt 1998): p, = Pg a 5 s /f,,, with U = 4 X 10 9 yrs. Stars are 
allowed to form in regions with temperature T < 2 x 10 4 K 
and gas density n > 0.1 crrT 3 . 5 The code also accounts for 
the stellar feedback on the surrounding gas, including injec- 
tion of energy and heavy elements (metals) via stellar winds, 
supernovae, and secular mass loss. The details of star forma- 
tion prescription and feedback are discussed in Kravtsov et al. 
(2005). Some potentially relevant physical processes, such as 
AGN bubbles, physical viscosity, magnetic field, and cosmic 
rays, are not included. 

The adaptive mesh refinement technique is used to achieve 
high spatial resolution in order to follow the galaxy formation 
and evolution self-consistently in these simulations. The peak 
spatial resolution is a 7 and 5/T 1 kpc, and the dark matter 

5 We have compared runs where star formation was allowed to proceed in 
regions different from our fiducial runs. We considered thresholds for star 
formation of n = 10, 1, 0.1, and 0.01 cirT 3 . We find that thresholds affect 
the properties of the ICM at small radii, r/r vlr < 0.1, but differences are 
negligible at the radii we consider in this study. 
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TABLE 1 

Simulated Cluster Sample of the CSF Run at z=0 



Name r 50Q Tx rel./unrel. 

(ft-'Mpc) (h- [ U) l3 M ) (keV) a (keV) b (I/O)' 

CL101 .... 1.160 8.17 90.81 7.44 8.72 8.67 8.86 000 

CL102 .... 0.978 4.82 54.47 5.63 5.63 5.83 5.86 000 

CL103.... 0.994 4.92 57.71 4.84 4.73 4.93 4.62 000 

CL104.... 0.976 5.15 53.88 6.61 7.69 7.73 7.73 111 

CL105.... 0.943 4.71 48.59 5.67 6.21 6.21 6.17 001 

CL106.... 0.842 3.17 34.65 4.54 4.34 4.35 4.30 000 

CL107.... 0.762 2.17 25.66 3.61 3.97 3.71 3.94 100 

CL3 0.711 1.91 20.90 3.37 3.65 3.60 3.61 111 

CL5 0.609 1.06 13.11 2.22 2.40 2.39 2.39 111 

CL6 0.661 1.38 16.82 2.88 3.38 3.38 3.57 000 

CL7 0.624 1.21 14.13 2.54 2.96 2.88 2.90 111 

CL9 0.522 0.73 8.23 1.58 1.53 1.60 1.57 000 

CLIO 0.487 0.43 6.72 1.58 1.93 1.90 1.91 111 

CL11 0.537 0.78 8.99 1.75 2.00 2.02 1.98 000 

CL14 0.509 0.62 7.69 1.64 1.85 1.84 1.83 111 

CL24 0.391 0.26 3.47 0.97 1.06 1.04 1.07 010 

is the average temperature measured directly from the 3D ICM distribution in the 
simulations. 

b Average temperatures measured in the shell of [0. 15,l]r5oo from the mock Chandra analysis 
of simulated clusters viewed along three orthogonal projection axes (x, y, z,from left to right). 
Note that the values of Tx quoted here are different from those in Table 1 of Nagai et al. (2007), 
where erroneous values were presented. 

c Classification of relaxed and unrelaxed clusters are indicated with and 1, respectively, for 
the three projections. 



particle mass in the region around the cluster was 9. 1 x 10 8 and 
2.7 x 10 8 h~ x M e for CL 101-107 and CL 3-24, respectively. 
To test the effects of galaxy formation, we also repeated each 
cluster simulation with only the standard gasdynamics for the 
baryonic component, without radiative cooling or star forma- 
tion. We will use labels "non-radiative" and "cooling+SF" 
(CSF) to refer to these two sets of runs, respectively. 

In this work, we also use mock Chandra X-ray images and 
spectra of the simulated clusters to derive total mass, gas mass 
and temperature profiles, as well as integrated cluster prop- 
erties, using the analysis procedures essentially identical to 
those used to analyze real Chandra observations, as described 
in Nagai et al. (2007). The average X-ray spectral temper- 
ature, T x , is obtained from a single-temperature fit to the 
spectrum integrated within r$oo, excluding the central region, 
r < 0.15r5oo- For each cluster, the mock data is created for 
three orthogonal projections along the x, y, and z coordinate 
axes. In § 5.3 we use quantities derived from the mock obser- 
vations to compare scaling relations exhibited by simulated 
clusters to observations. 

Our simulated sample includes 16 clusters at z — and 
their most massive progenitors at z = 0.6. The properties of 
simulated clusters at z = are given in Table 1 . The masses 
are reported at the radius rsm enclosing overdensities with re- 
spect to the critical density at the redshift of the output. This 
choice of the outer radius is mainly motivated by the fact that 
clusters are more relaxed within rsno compared to the outer 
regions (Evrard et al. 1996). We also use r2oo, nooo, and r25oo 
which are approximately 1.52, 0.71, and 0.44 times rsoo, re- 
spectively. Mean spectral temperatures are presented sepa- 
rately for the three orthogonal projections to show the varia- 
tion due to projection effects, substructure, etc. Note that the 
values of T x quoted here are different from those in Table 1 
of Nagai et al. (2007), where erroneous values were presented 
by mistake. In our analysis below we distinguish unrelaxed 
and relaxed clusters for a more consistent comparison with 
the observations. The classification is based on the overall 



morphology of the mock X-ray images, as discussed in Nagai 
et al. (2007). In Table 1 relaxed and unrelaxed clusters are 
indicated with and 1 for the three orthogonal projections (x, 
y, z from left to right). 

3. OBSERVATIONAL CLUSTER SAMPLE 

To test our simulation results against observations we use 
a set of accurate measurements of gas density, temperature, 
and total mass profiles for a sub-sample of 13 relaxed clus- 
ters at z ~ that was presented in Vikhlinin et al. (2005) 
and Vikhlinin et al. (2006). The clusters are selected on the 
basis of regular and relaxed morphology of their X-ray sur- 
face brightness images, although some of the systems show 
signs of AGN activity in their cores. Three of the low- 
T x clusters, including USGC S152 (T x = 0.69 keV), A262 
(T x = 1.89 keV), and RXJ1 159+5531 (T x = 1.80 keV), are 
excluded from the comparisons that involve measurements of 
M 5 oo or normalization with ^500 because their values are very 
uncertain due to insufficient spatial coverage. In § 5.2, we 
include USGC S152 and RXJ1 159+5531 for comparisons of 
the entropy scaling relations at 0.1r2oo, ^2500, and riooo, and 
A262 at the first two radii, but not at nooo- Since none of 
the measurements extends out to r2oo, we estimate r2oo us- 
ing ?"2oo = 1 -52r 5 oo, which provides a robust and accurate es- 
timate of r2oo for our CSF and non-radiative simulations as 
well as the XMM-Newton mass measurements Pointecouteau 
et al. (2005). The observations and analysis procedure used 
to extract ICM properties and profiles from the Chandra data 
are described in detail in Vikhlinin et al. (2005) and Vikhlinin 
(2006). 

In our previous study (Nagai et al. 2007), we used the mock 
Chandra images and spectra of the simulated clusters to as- 
sess the accuracy of the X-ray measurements of galaxy clus- 
ter properties. Our results show that the X-ray analysis of 
Vikhlinin et al. (2006) provides very accurate reconstruction 
of the 3D gas density and temperature profiles for relaxed 
clusters. Therefore, we directly compare the profiles derived 
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from Chandra analysis to the 3D profiles of simulated clus- 
ters. Note that masses and overdensity radii of the clusters in 
the Chandra sample were derived from the X-ray hydrostatic 
analysis. A bias in the estimated cluster mass may results in a 
slight underestimate of the estimated cluster virial radius r^J, 
by about a few percent for relaxed clusters (see also Nagai 
et al. 2007, for more details and discussions). We will show a 
such comparison in § 4. 

At the same time, our tests show that X-ray analysis can 
result in a a 15% underestimate in the hydrostatic estimates 
of total cluster mass. The bias is due to the non-thermal pres- 
sure support from the sub-sonic turbulent motions of the ICM 
gas, ubiquitous in cluster simulations (Evrard 1990; Norman 
& Bryan 1999; Nagai et al. 2003; Rasia et al. 2004; Kay et al. 
2004; Faltenbacher et al. 2005; Dolag et al. 2005; Rasia et al. 
2006; Nagai et al. 2007), but not included in observational 
hydrostatic mass estimates. In § 5.3, we present the com- 
parisons of the X-ray observable-mass relations of our simu- 
lated clusters to deep Chandra X-ray observations of nearby, 
relaxed clusters using both the true masses of clusters mea- 
sured in simulations and the masses estimated from hydro- 
static equilibrium analysis. We also correct for the differ- 
ences in the assumed cosmological parameters in simulations 
(fb = Qb/Qm = 0.1429 and h = 0.7) to those assumed in 
the observational analyses (fb = 0.175 and h = 0.72). We 
adopt fb = 0.175 and h = 0.7 throughout this work. Note that 
Qm = 0.3 and £2a = 0.7 are assumed in both analyses. 

In addition, we compare in § 5.2 the Chandra ICM entropy 
measurements with those based on the XMM-Newton observa- 
tions of 10 clusters (Pratt et al. 2006) and the ROSAT+ASCA 
data for 64 clusters (Ponman et al. 2003). In these compar- 
isons, we do not use measurements that involve extrapolation, 
so that we minimize biases arising from such procedure. 

4. EFFECTS OF GALAXY FORMATION ON THE ICM 
PROFILES 

In this section we investigate the effects of galaxy formation 
on the ICM properties by comparing simulations performed 
with and without the processes associated with the galaxy 
formation: gas cooling, star formation, stellar feedback, and 
metal enrichment. Similar comparisons have been done for 
a subset of 1 1 clusters in Kravtsov et al. (2005) and Nagai 
(2006). Here, we use the extended sample of 16 clusters and 
analyze the sub-sample of six relaxed clusters that are identi- 
fied as "relaxed" in all three orthogonal projections, indicated 
as 1 1 1 in the last column of the Table 1. In § 5, we compare 
the results of simulations to the Chandra X-ray observations 
of nearby relaxed clusters. 

Figure 1 shows the average radial profiles of the ICM in 
relaxed clusters at z = in the CSF and non-radiative runs. 
Clockwise from the top-left panel, we show the gas density, 
temperature, pressure, and entropy profiles. The mean pro- 
files are obtained by first normalizing the ICM profiles of each 
cluster at rsoo and then averaging over a sample of relaxed 
clusters. The shaded bands show lcr rms scatter around the 
mean profile of the CSF runs, and the mean and scatter of the 
profiles are computed for a logarithm of each thermodynamic 
quantity. We also examine systems with T x >2.5 keV and 
<2.5 keV separately to study the mass dependence of the ef- 
fects (and also the effects of cooling in the bremsstrahlung- 
and line emission-dominated regimes). In the bottom panel of 
each figure, we also show the fractional change of the ICM 
profiles in the CSF runs relative to the non-radiative runs. 

The temperature, entropy (defined as K = ksT/n 2 ^), and 



pressure profiles are normalized to the values computed for 
the given cluster mass using a simple self-similar model 
(Kaiser 1986; Voit 2005): 
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where M500 is a total cluster mass enclosed within ^00, 
E 2 (z) = Qm(1 + zf + Cl A for a flat universe with a cos- 
mological constant assumed in our simulations. Numerical 
coefficients in equations (2) and (3) follow from the defini- 
tions ^500 = k B T 500 /n^ 5 3 00 and P500 = « g ,500 k B T 500 , where 
« c ,500 = (j"/y"c) «g,500 = 500 fbPcm/(M e m p ), p crit is the critical 
density of the universe, fb = Qm/^ is the mean baryon frac- 
tion in the Universe, /i is the mean molecular weight, and // e 
is the mean molecular weight per free electrons. Note that we 
use \x = 0.59 and /j. e = 1.14 throughout this work. 

Figure 1 shows that including the gas cooling and star for- 
mation significantly modifies the ICM profiles throughout the 
cluster volume. The effect is larger in the inner region and for 
the systems with lower T x (or the cluster mass). Compared 
to the non-radiative runs, the gas density in the CSF runs is 
reduced by a 50 and 25% at 0.3 and 1.0 rsoo, because a frac- 
tion of gas is converted into stars. At small radii, we observe 
a trend with cluster mass; for example, the suppression of the 
gas density in the CSF runs at r = 0. 1 ^00 is ~ 50% and 30% 
for systems with Tx < 2.5 keV and > 2.5 keV, respectively. 
However, at r > OJrsoo (or r > 0.15r v ; r ), our simulations 
show very little systematic trend with T x , indicating that the 
clusters become self-similar in the outskirts even when the 
cooling and star formation are turned on. 

The ICM temperature profiles decline monotonically from 
0.05r5oo outwards. The shape of the temperature profiles are 
similar between the non-radiative and CSF runs, but there is a 
clear offset in their normalization. The temperature in the CSF 
runs is systematically higher by 10%-20% outside the core, 
indicating that the net effect of gas cooling and star formation 
is to increase the ICM temperature. The effects of gas cooling 
and star formation on the ICM temperature show a stronger 
dependence on cluster mass than gas density. 

The ICM entropy provides one of the most fundamental in- 
sights into physical processes determining the thermodynam- 
ics of the ICM, because it is expected to a be a conserved 
quantity, modified only by shock waves and "non-adiabatic" 
processes we are interested in (e.g., see Voit et al. 2002, 2005, 
and references therein). Figure 1 shows that the entropy pro- 
files in our non-radiative simulations scale self-similarly and 
are well-described by a power law K oc r 1 2 , at r > 0.3r5oo, in 
agreement with previous studies (Voit et al. 2005). Note, how- 
ever, that there is a systematic discrepancy between the pre- 
dictions of the Eulerian and SPH codes at small radii (Frenk 
etal. 1999; Ascasibaret al. 2003). However, since the primary 
focus of this paper is on the ICM properties outside the cluster 
core, we leave a detailed analysis of this entropy discrepancy 
for future work. 

Using the average K(r) profile from the non-radiative sim- 
ulations as a baseline, we study the effects of gas cooling and 
star formation on the ICM entropy. Compared to the non- 
radiative runs, the ICM entropy in the CSF runs is enhanced 
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Fig. 1 . — Radial profiles of the ICM in relaxed simulated clusters at z = 0. For each of the physical profiles the upper panels show the profiles, while the bottom 
panels show the corresponding fractional deviations of the profiles in the CSF simulation from the corresponding profiles in the non-radiative runs. The figure 
shows gas density (top-left), temperature (top-right), entropy (bottom-left), and pressure (bottom-right) profiles. Thick solid and dashed lines show the average 
profiles of the relaxed clusters in the CSF and non-radiative runs, respectively. The shaded band indicates the rms scatter around the mean profile for the CSF 
run. In addition, the dashed and dotted lines indicate the average profiles of systems with Tx > 2.5 and < 2.5 keV, respectively, in the CSF simulations. Note 
that the entropy profiles of the non-radiative runs outside 0.3r5oo are well-described by a power-law radial profile K oc r L2 , indicated by the dashed line in the 
bottom-left panel. 



in the entire radial range of interest, even at the virial radius 
r v i, =i 2 x r5oo- This is because cooling leads to the condensa- 
tion of the lowest entropy gas, which is replaced by the gas of 
higher entropy (Bryan 2000; Voit & Bryan 2001). The effect 
strongly depends on radius and is most pronounced in the in- 
ner regions with the largest effect of ~ 100% near r = 0.3r5oo. 
However, the entropy is enhanced by ~ 40% even at r — r^oo. 
The figure also shows that the magnitude of the effect in the 
inner regions depends on cluster mass, but is approximately 
the same for all masses at r > 0.5r5oo), indicating that cooling 
preserved self-similarity of the cluster outskirts. 

Finally, pressure profiles exhibit the most remarkable de- 
gree of self-similarity and low-level of cluster-to-cluster scat- 



ter. Notice that the average pressure profiles of low and high- 
Tx systems are nearly identical. This indicates that the self- 
similarity is best preserved for the quantities directly propor- 
tional to the ICM pressure or thermal energy, such as the inte- 
grated pressure F S z K ^g^mg (Nagai 2006) and Y x = M g T x 
(Kravtsov et al. 2006). Note, however, that inclusion of gas 
cooling and star formation modifies the overall shape and nor- 
malization of the pressure profiles and hence the Ysz and Yx 
parameters for the clusters of a fixed mass. In our simulations, 
the ICM pressure is suppressed by about 25% and 40% at 7500 
and r 2 5oo, respectively. 

5. COMPARISONS WITH OBSERVATIONS 
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Fig. 2. — Comparison of the ICM profiles in relaxed clusters at the present day (z a 0) in cosmological cluster simulations and the Chandra sample of Vikhlinin 
et al. (2006). The panels show the gas density (top-left), temperature (top-right), entropy (bottom-left), and pressure (bottom-right). Thick solid and dashed 
lines show the mean profiles in the CSF and non-radiative simulations, respectively, while the observed profiles are shown by the thin dotted, long-dashed and 
short-dashed lines for the systems with Tx > 5 keV, 2.5 < Tx < 5 keV, and Tx < 2.5 keV, respectively. Note that at r > 0.1 ("500 the profiles of the CSF 
simulations provide a better match to the observed profiles than the profiles in the non-radiative runs. 



In this section we present detailed comparisons of gas den- 
sity, temperature, entropy, and pressure profiles in the simu- 
lated clusters and the Chandra observations of low-z relaxed 
clusters (Vikhlinin et al. 2006). We also compare the scaling 
relations between Tx, M s , the X-ray "pressure" (Fx = M g Tx), 
and cluster total mass. 

5.1 . Profile comparison 

Figure 2 compares the ICM profiles in simulations with 
those observed. For the simulated clusters, we plot the mean 
ICM profiles averaged over the samples of relaxed clusters in 
both non-radiative and CSF runs. They are compared to the 
Chandra X-ray measurements of 1 1 nearby relaxed clusters. 
The observed clusters with Tx > 5 keV, 2.5 < Tx < 5 keV, 
and Tx > 5 keV are indicated with thin dotted, long-dashed, 



and short-dashed lines. 

The comparisons show that the ICM profiles in the CSF 
simulations agree quite well with observations outside the 
cores of clusters (r > f25oo)> while the non-radiative simu- 
lations predict overall shape and normalization of the ICM 
profiles inconsistent with observations at all radii. The obser- 
vations indicate then that a significant amount of cooling and 
condensation of gas out of hot ICM phase has occurred in real 
clusters. 

The ICM profiles in the inner regions, on the other hand, 
are not reproduced well in any of our simulations. The only 
exception is the n gas (r) profiles for Tx > 5 keV clusters, 
where the CSF runs produce results consistent with observa- 
tion down to r - 0.06 rsoo- However, even for these clus- 
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Fig. 3. — Correlation between the entropy K = kgTlnJ as a function 
of Mjog. The entropy scaling relations are measured at 0.1r2oo. ^soo, Hooo, 
and rsoo- We compare the relations in the CSF and non-radiative simulations, 
indicated with filled and open circles, and the dashed and dotted lines indicate 
the best-fit power law relations to these sets of simulations, respectively. Stars 
and triangles are observations by Chandra (Vikhlinin et al. 2006) and XMM- 
Newton (Pratt et al. 2006). The dashed lines indicate the best-fit power law 
relations to the CSF (upper line) and non-radiative (lower line) simulations, 
while the dotted and dot-dashed lines indicate the fits to the Chandra, and 
XMM-Newton measurements, respectively. 



ters the discrepancies between simulations and observations 
are evident in the T(r) and P(r) profiles at r < 0.3rgoo- The 
shape of the entropy profiles is different in the inner region 
of the simulated and observed clusters. For example, the en- 
tropy profiles of the observed clusters monotonically decrease 
with decreasing radius, while simulated entropy profiles flat- 
ten substantially at r < 0.2rsoo- 

Note that the ICM profiles of the T x < 2.5 keV Chan- 
dra clusters are systematically offset from the high-7x clus- 
ters in the inner regions and exhibit more pronounced cluster- 
to-cluster variations. The lowest-Tx system, MKW4 (Tx = 
1.64 keV), shows the most striking deviations from self- 
similarity. Our simulated clusters show similar trends, but 
the sample is too small to quantify the trends with Tx and the 
scatter. Note also that, if we use the estimated rsoo of simu- 
lated clusters to compare with data, the simulation curves in 
Figure 2 could slide to the left, bringing the characteristic val- 
ues plotted in Figures 3 and 4 down slightly (< 8%; see § 5.3 
for more discussions). 

5.2. Entropy scaling relations 

The scaling of the entropy with mass or Tx of clusters 
provides one of the most powerful diagnostics of the effects 
of galaxy formation on the ICM and deviations from self- 
similarity (Evrard & Henry 1991; Ponman et al. 1999; Voit 
& Bryan 2001; Pratt et al. 2006). For the self-similar cluster 
model, the entropy at a fixed overdensity radius is expected to 
scale linearly with Tx and with mass as oc M 2/3 (cf. eq. 1). 
In the discussion above we showed that entropy profiles of 
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Fig. 4. — Correlation between the entropy K = ksT/n^ as a function 
of Tx- The entropy scaling relations are measured at O.lr^oo, ^soo, nooo, 
and T5oo- Solid circles indicate the CSF simulations, while stars, triangles, 
and crosses are Chandra (Vikhlinin et al. 2006), XMM-Newton (Pratt et al. 
2006), and ROSAT+ASCA (Ponman et al. 2003) measurements. The dashed, 
dotted, and dot-dashed lines indicate the best-fit power law relations to the 
CSF simulations, Chandra, and XMM-Newton measurements, respectively. 



both observed and simulated clusters become approximately 
self-similar outside the cluster cores. In this section we ex- 
plicitly consider the scaling of entropy with cluster mass and 
temperature at four different radii and compare the results of 
numerical simulations with X-ray measurements obtained us- 
ing Chandra (Vikhlinin et al. 2006). We also compare the 
Chandra measurements with the XMM-Newton (Pratt et al. 
2006) and ROSAT+ASCA (Ponman et al. 2003) results. 

The entropy levels measured in the simulated clusters at 
0. 1 r2oo, ^2500; '"iooo, and rsoo are shown as a function of M500 in 
Figure 3 (such a relation was first studied observationally by 
Pratt et al. 2006). The best-fit power law approximations of 
these data, E(z) 2/3 K = C(M 500 /2 x lO 14 /r'M ) a , are given 
in Table 2 (the E(z) term corrects for evolution in the self- 
similar model, which needs to be applied to clusters observed 
at z + 0). Clearly, in the simulated clusters, inclusion of gas 
cooling and star formation increases the entropy level, and 
the magnitude of the effect is larger in the inner region. The 
changes in the normalization for M — 2 x 10 14 /i _1 M G clus- 
ters is a factor of 1.71, 1.65, 1.41, and 1.34 at r = 0.1 r 2 oo, 
7"2500> nooo, and rsoo, respectively. The effects of cooling and 
star formation on the ICM entropy are thus stronger at small 
radii. The slopes are consistent with the prediction of the 
self-similar model (a = 2/3) within 1 cr at all radii, in both 
CSF and non-radiative simulations. However, there are in- 
dications that the slopes are slightly shallower than the self- 
similar value, and the slope in the non-radiative run is some- 
what steeper than that in the CSF at r < rsoo- A larger sample 
of simulated clusters is needed to determine whether these 
differences are real. 

It is easier to compare the entropy levels in the simulated 
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TABLE 2 

Best-fit parameters for the entropy-mass relation, 
E(z) 2 ' 3 K = C(M 500 /2 X 10 14 /T 1 M o ) a at z = 0. 



radius 


quantity 


cooling+SF 


non-radiative 


Chandra 


XMM-Newton 


0.1r 2 oo 


loj 


oc 


2.423 ± 0.030 
0.647 ± 0.077 


2.189 ±0.025 
0.774 ± 0.062 


2.391 ±0.020 
0.335 ± 0.046 


2.286 ± 0.033 
0.304 ± 0.080 


("2500 


loj 


X 


2.734 ± 0.017 
0.593 ± 0.046 


2.516 ±0.012 
0.638 ±0.031 


2.805 ±0.011 
0.492 ± 0.038 


2.807 ± 0.014 
0.407 ± 0.044 


nooo 


lot 


IwC 
oc 


2.890 ± 0.016 
0.570 ± 0.043 


2.741 ±0.017 
0.596 ± 0.039 


2.984 ±0.011 
0.569 ±0.018 


3.094 ±0.013 
0.390 ± 0.046 


^500 


lot 


ho c 
oc 


3.063 ± 0.013 
0.605 ± 0.036 


2.936 ±0.013 
0.601 ± 0.037 


3.100 ±0.013 
0.598 ± 0.028 





TABLE 3 

Best-fit parameters for the entropy-temperature relation 
E(zf /3 K = C(T x /5 keV) a at z = 0. 



radius 


quantity 


cooling+SF 


Chandra 


XMM-Newton 


0.1r 2 oo 


logi C 

OC 


2.593 ± 0.042 
0.954 ± 0.077 


2.437 ± 0.019 
0.502 ± 0.084 


2.320 ± 0.035 
0.506 ±0.131 


^2500 


Jogio C 

(X 


2.905 ± 0.009 
0.958 ±0.031 


2.876 ±0.013 
0.769 ± 0.054 


2.851 ±0.021 
0.657 ± 0.089 


nooo 


log io C 

X 


3.064 ± 0.006 
0.922 ± 0.026 


3.067 ± 0.015 
0.887 ± 0.043 


3.135 ±0.026 
0.616 ±0.106 


'"soo 


log 10 C 


3.232 ±0.018 
0.946 ± 0.060 


3.186 ±0.018 
0.930 ± 0.063 





and observed clusters via the K - T x correlation. The results 
for our CSF runs are r = 0.1 r2oo, ^soo, Hooo. and r 500 are 
shown in Figure 4. Table 3 lists the best-fit parameters of the 
power approximations, E(z) A ^K = C(Tx/5 keV) a . As in the 
case of the K - M relation, the power-law slopes in the CSF 
runs are slightly shallower than, but consistent with, the self- 
similar expectation (a = 1) at all radii considered. Note also 
that the K - Tx relations exhibit remarkably tight relations at 
r > >"2500 for all clusters. 

We also show in Figure 3 and 4 the K - M and K - 77 scal- 
ings derived from several sets of X-ray cluster observables. It 
is most straightforward to compare our simulations with the 
Chandra results of Vikhlinin et al. (2006) (shown by stars in 
Figures 3 and 4), because we explicitly tested their data anal- 
ysis procedures (Nagai et al. 2007) and because the Chandra 
results for many clusters extend to r^oQ. First, we note the 
entropy normalizations in the Chandra clusters show a good 
overall agreement with the CSF runs at all radii. The results 
for non-radiative runs are strongly inconsistent with the data 
(see also Ponman et al. 2003; Pratt et al. 2006). 

At small radii, r25oo and 0.1r2oo, the observed clusters show 
significantly shallower slopes than expected in the self-similar 
model and seen in the simulations. For example, the slope of 
the K - 77 relation at r = 0.1r 2 oo is «csf = 0.95 + 0.08 for 
the simulated clusters, and achandra = 0.50 + 0.08 for Chandra 
sample. These results are in line with the disagreement be- 
tween the K(r) profiles of the simulated and observed clusters 
at small radii, discussed in § 5.1. 

The agreement, however, improves quickly as we go to 
larger radii. At r2soo, the slopes of the K - 77 relation are 



«csf = 0.96 + 0.03 for simulated clusters and occhandra 
0.77 ± 0.05 for the Chandra sample. At larger radii, nooo and 
rsoo, the Chandra-observed relations are fully consistent with 
the CSF simulations both in terms of slope and normalization; 
the slopes are also very close to the self-similar expectations, 
olk-m = 2/3 and ct K -T = F These results confirm the gen- 
eral conclusion of § 5.1 that although effects of cooling on the 
entropy normalization are significant within radii as large as 
rsoo, the scaling of the thermodynamic properties of the ICM 
become close to the self-similar expectation outside the inner 
cluster region. 

Also shown in Figure 3 and 4 are the entropy scaling re- 
lations derived from two more X-ray data sets, the XMM- 
Newton sample of Pratt et al. (2006) [triangles], and the 
ASCA+ROSAT sample of Ponman et al. (2003) [crosses at 
r = 0.1 f2oo]- At small radii, where the X-ray measurements 
are most straightforward, there is a good agreement between 
all observed relation. In particular, the XMM-Newton and 
Chandra relations for r = 0. 1 r2oo and r25oo are nearly identi- 
cal. A small offset of the Ponman et al. (2003) data points can 
be explained by the slightly different definitions of r2oo used 
in these works 6 . 

At r = riooo, the entropy normalizations for the most mas- 
sive clusters are in agreement for the Chandra and XMM- 
Newton samples but there is some tension in the values of 
slope. The XMM-Newton results indicate nearly the same 
slopes at r — nooo and smaller radii, all significantly flat- 
ter than the self-similar prediction: (Xk-t = 0.51 ± 0.13, 
0.66 ± 0.09, and 0.62 ± 0.1 1 for r = 0.1 r2oo, ^"2500, and nooo, 
respectively. The Chandra results clearly indicate a signifi- 
cantly steeper slope at this radius, 0.89 ± 0.04; a steep slope, 
0.93 + 0.06 is also observed at r — r^m- The statistical sig- 
nificance of the difference between the Chandra and XMM- 
Newton slopes is 3.4cr for the K - M and 2.4 cr for the K -T 
relations. The best-fit values of a indicate qualitatively differ- 
ent cluster properties. While the XMM-Newton results suggest 
that the departures from self-similar scalings are of the simi- 
lar amplitude at all radii, the Chandra measurements clearly 
point in the direction that the effect is confined to the very 
central regions. 

A detailed comparison of the XMM-Newton and Chandra 
data analyses is beyond the scope of this work. We only point 
out two effects that may contribute to the difference in the en- 

6 The quantity r2oo was defined in Ponman et al. (2003) through the Evrard 
et al. (1996) scaling with Tx', through the NFW fit to the data at smaller 
radii in Pratt et al. (2006); and through hydrostatic estimates of r^oo for the 
Chandra clusters. 



EFFECTS OF GALAXY FORMATION ON ICM 



9 



tropy scalings. First, the Chandra temperature profiles show 
a systematic decline at large radii (by a factor of ~ 1.7 at 
r = 7"5oo relative to peak values near ~ 0.2 rson), the XMM- 
Newton temperature profiles in the Pratt et al. (2006) sam- 
ple are much flatter. This systematic difference is discussed 
in Vikhlinin et al. (2005). Second, the gas densities in the 
Chandra analysis were derived using a model that allows for 
steepening of the p(r) profile at large radii. The XMM-Newton 
data were fit with the /?-type models that do not allow for such 
steepening (Pointecouteau et al. 2005). This leads to some- 
what different gas density profiles at large radii (see, e.g., Ap- 
pendix A2 in Vikhlinin et al. 2006). 

5.3. Relations between Total Mass and X-ray Observables 

We present comparisons of the X-ray observable-mass rela- 
tions of the CSF simulations and Chandra X-ray observations 
of nearby relaxed clusters in Figure 5. Following Kravtsov 
et al. (2006), we consider three X-ray proxies for the cluster 
mass: — the spectral temperature (7x), the gas mass (M g ), 
and the X-ray pressure (Yx = TxM g ). These X-ray mass 
proxies are derived from mock Chandra images of the sim- 
ulated clusters and analyzing them using a model and pro- 
cedure essentially identical to those used in real data anal- 
ysis. Note that the mean temperatures were estimated from 
a single-temperature to the spectrum integrated in the radial 
range [0.15 - l]r 5 oo (i.e., excluding emission from cluster 
core). 

In the upper panels of Figure 5, we compare the scal- 
ing relations of simulated clusters for the true cluster mass, 
M^™(< ^qq), measured in simulations to the relations from 
the Chandra X-ray cluster observations. We also plot the 
best-fit M-Tx relation from the XMM-Newton measurements 
(Arnaud et al. 2005) for comparison. Results of power law fits 
to these relations for different subsets of the clusters are pre- 
sented in Table 2 of Kravtsov et al. (2006). 7 The comparisons 
show that the normalizations of the scaling relations involving 
true M500 for our simulated sample are systematically high by 
~ 10 - 20% compared to the observed relations. We note that 
this level of agreement is considerably better than agreement 
between simulations and data as recently as several years ago 
(e.g., Pierpaoli et al. 2001). 

The remaining bias could arise from the assumption of the 
hydrostatic equilibrium, which is a key assumption that en- 
ables measurements of gravitationally bound mass of clusters 
from the X-ray and Sunyaev-Zel'dovich (SZ) effect data. To 
illustrate this, we compare the scaling relations based on the 
estimated hydrostatic mass (M^, E (< r 5o )) derived from the 
mock Chandra analysis of simulated clusters to the Chandra 
measurements in the lower panels of Figure 5. Note that we 
account for additional biases in the estimated cluster masses 
arising from a bias in the estimation of a cluster virial radius 
by measuring the cluster mass within the r|^, estimated from 
the hydrostatic analysis (see also Nagai et al. 2007, for more 
details and discussions). Similarly, the gas mass (M gi 5oo) is 
also computed at r^™ and r?5L in the upper and lower pan- 
els. 8 In Table 4, we summarize results for the relaxed clusters 
at z=0, relevant for comparison with observations considered 
here. These analyses show that the simulation results are in 
much better agreement with observations when using the hy- 



TABLE 4 

Best-fit parameters and scatter in the mass vs. proxy relations, 
m5q0 = cx 01 , for relaxed clusters at 2 = and 0.6. 



relation 8 


quantity b 


M 5 s ™(True) 


M 5 S M(HSE) 


M ^(Chandra) 


M500 


-T X 


log 10 C 
scatter 

lOglo C ss 


14.21 ±0.010 

1 < 1 < _i_ n nco 
l.JlJ ± U.Ujz 

0.136 
14.21 ± 0.009 


14.10 ±0.008 

l.OZi ± U.UZ/ 

0.117 
14.10 ±0.009 


14.11 ±0.035 
1 A qq _i_ n noi 

1 .4sy ± u.uyj 
14.10 ±0.014 


M500 


-M g 


logioC 

a. 
scatter 
log 10 C ss 


14.12 ±0.008 
0.894 ± 0.023 

0.114 
14.12 ±0.010 


14.05 ±0.011 
0.932 ± 0.033 

0.153 
14.05 ±0.011 


14.14 ±0.044 
0.811 ±0.067 

14.07 ± 0.022 


M500 


-¥ X 


logioC 

a. 
scatter 
logio C ss 


14.06 ± 0.004 
0.568 ± 0.006 

0.053 
14.05 ± 0.005 


13.97 ± 0.007 
0.596 ±0.010 

0.087 
13.97 ±0.006 


14.04 ± 0.047 
0.526 ± 0.038 

13.98 ±0.017 



"Power law fits were performed for relaxed clusters of our cluster sample at 
z = and 0.6. In addition to the fits in which both normalization and slope of 
the power-law relations were fit simultaneously, we provide the best fit normal- 
izations, C ss , for each relation when fit with the slopes fixed to their self-similar 
values: 1.5, 1.0, and 0.6 for the M s00 -Tx, M s00 -M g , and M 500 -Yx relations, 
respectively. 

b For each observable X (= Tx, M g , Yx), we fit power law relation of the form 
M s00 = CiX/Xo) 01 , with X = 3.0 keV, 2 X 10 13 M Q , and 4 X 10 13 keV M Q , for 
Tx, M g , Yx, respectively. Note that M500 is in units of (/r'M ). 



7 Note that M500 is in units of h 
Kravtsov et al. (2006). 



Mo in this work, while it is Mm 



Note that Tx is computed within in both panels; however, correcting 
for the bias in r 500 has a negligible (< 1%) effect on the Tx estimate. 



drostatic mass. The systematic offset in normalizations could 
thus be due to the bias of total hydrostatic mass estimate due 
to turbulent motions of the ICM. 

In terms of scatter, the M 5u o-7x relation exhibits the largest 
scatter of ~ 20% in M 500 around the mean relation, most of 
which is due to unrelaxed clusters. The unrelaxed clusters 
also have temperatures biased low for a given mass, likely be- 
cause the mass of the system has already increased but only a 
fraction of the kinetic energy of merging systems is converted 
into the thermal energy of gas, due to incomplete relaxation 
during mergers (Mathiesen & Evrard 2001). Unfortunately, 
we cannot compare the scatter directly to the Chandra results 
because for real clusters, the scatter is dominated by the mea- 
surement uncertainties and the intrinsic scatter (see discussion 
in Vikhlinin et al. 2006). The slope and redshift evolution 
of the M500 - T x relations are quite close to the simple self- 
similar expectation. The M500 - M g relation has a somewhat 
smaller scatter (a 11%) around the best-fit power-law relation 
than the M500 - Tx, but its slope is significantly different from 
the self-similar prediction for the M 5 no - M g relation due to 
the trend of gas fraction with cluster mass present for both the 
simulated clusters in our sample (see Kravtsov et al. 2005) 
and for the observed clusters (Vikhlinin et al. 2006). In all of 
the scaling relations considered here, the use of M^. E , instead 
of M^, modifies the scatter by a few percent for the relaxed 
clusters. 

The Msqq-Yx relation shows the scatter of only « 7%, mak- 
ing it by far the tightest of all the scaling relations. Note that 
this value of scatter includes clusters at both low and high- 
redshifts and both relaxed and unrelaxed systems. The tight- 
ness of the M500 - ix relation and simple evolution are due 
to a fortunate cancellation of opposite trends in gas mass and 
temperature (see Kravtsov et al. 2006). The slope and red- 
shift evolution of normalization for the M500 - Yx relations 
are well-described by the simple self-similar model. 

6. DISCUSSION AND CONCLUSIONS 
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Fig. 5. — Comparisons of the X-ray observable-mass relations in simulations and observations. From left to right: Correlations between the total mass, M500 
and X-ray spectral temperature (Tx), gas mass (A/ e ), and X-ray pressure (Yx 



measured in simulations (upper panels) and the hydrostatic mass Ml 



HSE . 
500 



TxM g ). Relations are shown for the true 3D cluster mass M500 = M(< r"^) as 



M (< rSL) derived from mock Chandra analysis (lower panels). Separate symbols 



indicate relaxed and unrelaxed simulated clusters, and also z = and 0.6 samples. The figures include points corresponding to three projections of each cluster. 
The dot-dashed lines are the power law relation with the self-similar slope fit for the sample of relaxed clusters. The dotted lines indicate the rms scatter around 
the mean relation. The data points with error bars are Chandra measurements of nearby relaxed clusters. The dashed line is the best-fit M-Tx relation from the 
XMM-Newton measurements. 



We presented analyses of the simulated cluster sample of 
16 clusters spanning a representative mass range (5 x 10 13 - 
2 x 10 15 hT x M ) and modeled using the shock-capturing Eu- 
lerian adaptive mesh refinement iV-body+gasdynamics ART 
code in the ACDM cosmology. These simulations achieve 
high spatial resolution and include various physical processes 
of galaxy formation, including radiative cooling, star forma- 
tion and other processes accompanying galaxy formation. We 
study the effects of gas cooling and star formation on the ICM 
properties by comparing two sets of simulations performed 
with and without these processes included. The results of sim- 
ulations with dissipation are compared to recent Chandra X- 
ray observations of nearby, relaxed clusters (Vikhlinin et al. 
2005, 2006). 

We show that gas cooling and star formation modify both 
the normalization and the shape of the gas density, tempera- 
ture, entropy, and pressure profiles. As the lowest-entropy gas 
cools and condenses out of the hot phase in the cluster pro- 
genitors, the gas density in their inner regions is lowered and 
entropy is increased (Bryan 2000; Voit & Bryan 2001). The 
effects have strong radial dependence and are the strongest in 
the inner regions r < O.lrsoo- At these inner radii simulation 
profiles do not match the observations. On the other hand, at 
r > O.lrsoo the profiles in the CSF simulations and observa- 
tions agree quite well, while profiles in the non-radiative runs 



disagree with observations at all radii within r^oo- 

In particular, the simulations with cooling can explain the 
observed high levels of entropy in observed clusters compared 
to the non-radiative expectation pointed out previously (Pon- 
man et al. 2003; Pratt et al. 2006). At r > risoo the cluster pro- 
files are approximately self-similar within current statistical 
error bars, while there is an indication that their best-fit slopes 
are slightly shallower than the self-similar value. Moreover, 
the slope and normalization of the entropy scaling relations 
in the simulated clusters are in good agreement with Chandra 
observations at riooo and r^oo, while the observed relations ex- 
hibits deviations from the self-similarity at r < r25oo- Note 
also that the results of Chandra and XMM-Newton measure- 
ments agree quite well within r < ^2500, but the significant dis- 
agreement is seen at fiooo- The statistical significance of the 
difference between the Chandra and XMM-Newton slopes is 
3Ao- for the K-M and 2 Act for the K-T relations. The differ- 
ence is critical for theoretical interpretation and implications, 
as the XMM-Newton scaling was billed as a major evidence 
for deviations from self-similar scalings at large radii (indeed 
at nooo, where XMM-Newton measurements is several sigma 
away from the slope of 1.0). Chandra results show that devi- 
ations at larger radii, if they exist, are much smaller. Despite 
the deviations from self-similarity in the ICM entropies, we 
show that pressure profiles in particular, show a remarkable 
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degree of self-similarity and exhibit very small scatter. 

We also present comparisons of scaling relations of cluster 
X-ray observables with total cluster mass in the simulations 
with cooling and recent deep Chandra observations. Specif- 
ically, we compare correlation of the spectral X-ray temper- 
ature, ICM gas mass (M g ), and the X-ray equivalent of in- 
tegrated pressure {Yx = M g Tx)- In these comparisons X- 
ray observables for the simulated clusters are derived from 
mock Chandra analysis using procedure essentially identical 
to those used in real data analysis. 

The slope and normalization of the M - T x and M - Y x 
relations in simulations and Chandra observations are in 
good agreement, and they are consistent with the simple self- 
similar expectation. In terms of scatter, the M500 - Yx relation 
shows scatter of only « 7%, making it by far the tightest of all 
the scaling relations. Note that this value of scatter includes 
clusters at both low and high redshifts and both relaxed and 
unrelaxed systems. The M500 - 7x relation, on the other hand, 
exhibits the largest scatter of ~ 20% scatter in M500 around 
the mean relation, most of which is due to unrelaxed clusters. 
The unrelaxed clusters also have temperatures biased low for 
a given mass, likely because the mass of the system has al- 
ready increased but only a fraction of the kinetic energy of 
merging systems is converted into the thermal energy of gas, 
due to incomplete relaxation (Mathiesen & Evrard 2001) dur- 
ing mergers. 

Moreover, these comparisons show that the normalizations 
of the scaling relations of relaxed clusters in simulations and 
observations agree at a level of about « 10% - 20%. This is 
a considerable improvement, given that significant disagree- 
ment existed just several years ago (see Finoguenov et al. 
2001; Seljak 2002; Pierpaoli et al. 2003). The residual sys- 
tematic offset in the normalization is likely caused by non- 
thermal pressure support from subsonic turbulent gas motions 
(Evrard et al. 1996; Rasia et al. 2004, 2006; Faltenbacher 
et al. 2005; Dolag et al. 2005, E. Lau et al. 2007, in prepara- 
tion). This contribution is approximately independent of clus- 
ter mass (Vazza et al. 2006, E. Lau et al. 2007, in preparation) 
and is not accounted for in X-ray hydrostatic mass estimates. 
For example, when we repeat the comparison of scaling rela- 
tions using hydrostatic mass estimates for the observed clus- 
ters, we find excellent agreement in normalizations, demon- 
strating explicitly that there is a systematic « 10% - 20% 
offset between hydrostatic mass estimate and the true mass in 
simulated clusters. 

Part of the non-thermal pressure support may also be con- 
tributed by cosmic -rays and magnetic fields. In practice, it 
may be difficult to distinguish between different sources of 
non-thermal pressure. A possible test is their radial depen- 
dence. Turbulent motions, for example, are in general smaller 
at smaller radii and the turbulent pressure gradient is corre- 
spondingly smaller. In the case of turbulent motions, we can 
therefore expect that the bias in the total mass estimate should 
decrease at smaller radii. This may not be the case for some 
other sources of non-thermal pressure, although recent mod- 
els of cosmic-rays contribution to the total pressure show a 
qualitatively similar radial dependence as the turbulent pres- 
sure (Pfrommer et al. 2007). 

The much improved agreement of the scaling relations and, 
especially, normalization and shape of the gas profiles be- 
tween simulations with cooling and star formation and ob- 
servations show that inclusion of galaxy formation in cluster 
simulations results in more realistic modeling of the hot ICM. 
This may indicate that gravitational dynamics and the basic 



cooling of the hot gas accompanying galaxy formation are 
the dominant processes determining thermodynamics of the 
ICM outside the cluster cores, while other processes, such as 
feedback, thermal conduction, viscosity, and cosmic -rays, are 
playing only a minor role for a large fraction of the ICM mass. 

Note, however, that the agreement between our simulations 
and observations is achieved by condensation of a significant 
fraction of hot gas into cold dense phase, which is subse- 
quently converted into stars. Thus, simulated clusters have 
« 40% of their baryons within r 5 oo in stellar form at z = 0, 
while the rest of the baryons are in the hot phase. Although 
the low hot gas mass fractions (a; 60% - 70% of the universal 
value) are consistent with observations (Vikhlinin et al. 2006; 
McCarthy et al. 2007), the high stellar fractions are not. Note, 
however, that the reduced stellar fraction with more efficient 
stellar or AGN energy feedback generally results in the pro- 
files that are in between those of the CSF and non-radiative 
runs and disrupts the good agreement between models and 
data. 

Observational estimates of the stellar mass fractions in 
groups and clusters range from a 5% - 10% (Eke et al. 2005) 
to* 15% - 20% (Lin et al. 2003; Gonzalez et al. 2007) of 
the universal baryon fraction or at least a factor of 2-3 lower 
than the fractions found in our simulations. This is a well- 
known discrepancy often referred to as the "overcooling prob- 
lem". Our results show that the X-ray and optical observa- 
tions appear to give seemingly contradicting constraints. The 
low observed stellar fractions imply existence of an efficient 
mechanism suppressing star formation in real clusters, while 
observed properties of hot ICM are not consistent with small 
amounts of cooling (i.e., predictions close to the non-radiative 
limit). At present it is not clear how these two observational 
constraints can be reconciled. 

Our tests indicate that profiles and average quantities (i.e., 
gas fractions) derived from analyses of modern X-ray data 
are robust and do not suffer any obvious biases (Nagai et al. 
2007). On the other hand, there are certain systematic un- 
certainties in estimates of stellar mass from optical observa- 
tions related both to possible low-surface stellar component 
missed in shallow observations (e.g., Gonzalez et al. 2005, 
2007; Lauer et al. 2007; Seigar et al. 2007) and to the un- 
certainties in the stellar population modeling of the observed 
photometry. It is unlikely, however, that any single uncer- 
tainty is large enough to account for the entire factor of 2-3 
discrepancy between stellar fractions in simulations and ob- 
servations. We note also that stellar fraction predicted by sim- 
ulations depends on the implementation of the feedback pro- 
cesses in simulations (e.g., Borgani et al. 2006). However, the 
current implementations of the feedback processes efficient in 
significantly suppressing stellar fraction are essentially ad hoc 
and it is uncertain whether the feedback is actually as efficient 
in practice. 

The progress in our understanding of these issues should 
come from detailed convergence studies and comparisons of 
simulation results done using different numerical codes, fur- 
ther comparisons of simulations with deep X-ray observa- 
tions, deeper observations, and thorough analysis of uncer- 
tainties in the optical estimates of cluster stellar masses. 
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TABLE A 1 

Best-fit parameters of the pressure profile. 



Runs Sample P c 5m (a,/3,-y) 



Chandra Obs. Relaxed 3.3 1.8 (1.3,4.3,0.7) 

Cooling+SF Relaxed 3.3 1.8 (1.3,4.3,1.1) 

Cooling+SF Unrelaxed 2.0 1.5 (1.4,4.3,1.0) 

Non-radiative Relaxed 38.0 3.0 (1.1,4.3,0.3) 

Non-radiative Unrelaxed 3.0 1.5 (1.4,4.3,0.9) 
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APPENDIX 
ANALYTIC PRESSURE MODEL 

The Sunyaev-Zel'dovich (SZ) effect is a direct probe of thermal energy content of the Universe and provides a unique and 
powerful probe of the structure formation and cosmology in the near future. The SZ observations are probing the integrated 
pressure of the ICM. Accurate analytic parameterizations of the ICM pressure profiles can therefore be useful for developing 
efficient cluster detection algorithm for upcoming SZ cluster surveys, analysis and interpretation of SZ effect observations, as 
well as theoretical modeling of cluster ICM. The fact that the self-similarity is best preserved for the pressure profiles and their 
low cluster-to-cluster scatter (see § 4) provides further motivation for the use of accurate pressure profile parameterizations. 

Here we present a simple analytic model of the pressure profile that closely matches the observed profiles of the Chandra X-ray 
clusters and results of numerical simulations in their outskirts. Since the gas pressure distribution is primarily determined by the 
gravitationally dominant dark matter component, we parameterize the pressure profile using the generalized NFW model, 

™ - Po (Al) 

P500 x^(l+x«)^>/ a 

where x = r/r s , r s = rsoo/csoo, P500 is given by equation 3, and (oc.,/3, y) are the slopes at r ~ r s , r » r s , and r <K r s , respectively. 
We find that a model with Po=3.3, C500 = ^5oo/^=l-8, and(a,j6,y) = (1.3, 4. 3, 0.7) provides a good description of the pressure 
profiles of the high-rx Chandra clusters within the observed range (r < rsoo) as well as the profiles of simulated clusters in 
0.5 < r/rsoo < 2.0. Figure Al shows generalized NFW fits to the pressure profiles of relaxed clusters in simulations and Chandra 
observations. In the outskirts, we set the slope to be y — 4.3, which is the average best-fit values for both non-radiative and CSF 
simulations. For the Chandra clusters with Tx > 5 keV, the inner slopes of the pressure profile are « 0.7. The inner slopes appear 
to be shallower for the lower- Tx systems, but they also show much larger cluster-to-cluster variation. Table Al summarizes the 
best-fit model parameters for observed and simulated clusters. For the relaxed systems, the same set of parameters with a different 
inner slope y = 1.1 produces the pressure profile of the CSF run, while that of the non-radiative run requires a very different set 
of parameters. The pressure profiles of the unrelaxed systems are generally less concentrated (smaller value of C500) with slightly 
different inner and outer slopes. 
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Fig. Al. — Generalized NFW fits to the pressure profiles of relaxed clusters in simulations and Chandra X-ray observations. Open circles and squares show the 
mean profiles in the CSF and non-radiative simulations, respectively. Thin dotted shows Chandra X-ray clusters with Tx > 5 keV. Thick lines show the best-fit 
generalized NFW model to simulation and observed profiles (see Table Al for the best-fit parameters). 
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